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We expand on previous work that treats relaxation physics of low-lying excited states in ideal, 
single electron, silicon quantum dots in the context of quantum computing. These states are of 
three types: orbital, valley, and spin. The relaxation times depend sensitively on system parameters 
such as the dot size and the external magnetic field. Generally, however, orbital relaxation times 
are short in strained silicon (10 -7 to 10 -12 s), spin relaxation times are long, (10 -6 to ^> Is), while 
valley relaxation times are expected to lie in between. The focus is on relaxation due to emission 
or absorption of phonons, but for spin relaxation we also consider competing mechanisms such as 
charge noise. Where appropriate, comparison is made to reference systems such as quantum dots in 
III-V materials and silicon donor states. The phonon bottleneck effect is shown to be rather small 
in the silicon dots of interest. We compare the theoretical predictions to some recent spin relaxation 
experiments and comment on the possible effects of non-ideal dots. 

I. INTRODUCTION 

The spin of an electron in silicon may act as an information carrier in future information technologies, from quantum 
computers to spintronics. For quantum information applications, the spin of cold, localized electrons in silicon can 
make a qubit with a low memory-error-rate due to the purifiability of the spin environment (a spin-0 nuclear isotope 
is available) and silicon's inherently weak spin-orbit interaction, which isolates information stored in the electron 
spin from charge movement and other noise. Quantum dots, in addition, allow for ready tunability and alignment of 
the confined electron (for physical transport, computation via qubit-qubit coupling, readout, initialization) the 
potential for a fabrication route with present-day lithographic techniques, and the enabling of a fast, DC-controlled 
two-qubit gate based on Heisenberg exchange |2J. The goal of designing and constructing quantum computers based 
on quantum dots requires characterization of all their physically relevant properties. For this, it is necessary to have a 
full toolbox of experimental diagnostics - in this case, one electron excited- state lifetime measurements as a function 
of external parameters such as temperature and magnetic field, and to be able to interpret these measurements in the 
light of theory. 

In this paper we present calculations of the dominant spin relaxation processes in ideal silicon quantum dot spin 
qubits |3J along with calculations for orbital and valley relaxation. Our particular emphasis is decay due to phonon 
coupling, since this mechanism is believed to be responsible for T\ across wide parameter regimes relevant to quantum 
information applications. We will indicate in detail how it can be distinguished experimentally from other mechanisms. 

A full, experimentally verified theory of the energy relaxation processes of the excited electronic states in silicon 
quantum dots is important for several reasons. First, it accomplishes a major step on the experimental path to 
determining the quantum coherence times of isolated spins in silicon (a preeminent goal in the verification of a qubit). 
Second, it helps corroborate our understanding of the material system and better characterizes the device under 
scrutiny; an example would be transport spectroscopy of nearby energy levels and their line widths. Third, it provides 
necessary parameters needed for the design of future experiments, systems, and architectures. Indeed, if the dominant 
spin qubit relaxation mechanisms are as we predict, we can not only validate the T\ lifetimes of silicon qubits, but also 
retrieve the relative magnitudes of the dominant spin-orbit coupling contributions inherent in the device (important for 
both silicon quantum computing and silicon spintronics applications), as well as the nature and energies of the various 
states above the two spin-qubit states. Finally, the lifetimes of excited orbital states are relevant to optical pumping 
schemes, many-phonon decoherence calculations, transport spectroscopy, beyond single-spin qubit implementations 
[4j[5], and other areas of quantum control. 

It has long been known that localized spins in silicon can have exceedingly long lifetimes, even at relatively high 
temperatures (> 2K) |6J. Theoretical predictions of spin decoherence times are notoriously difficult as many mecha- 
nisms can relax a spin, even in the more robust case of direct energy relaxation, or T\ processes, where a quantum 
of energy is lost to the environment. In silicon, for example, energy relaxation processes may depend on the many- 
valley nature of the conduction band electrons; neglecting this effect leads to incorrect predictions (to many orders 
of magnitude) [7J[8]. The key realization came in 1960 from Roth [10] and, independently, Hasegawa [llj - that 
spin mixing to the (Is- like) valley manifold states explains the "fast" relaxation observed for donor electrons. Soon 
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after, Feher, Gere, Wilson et al HI] and Castner [131 E] thoroughly fleshed out the experiments and theories of 
donor state lifetimes in silicon. Castner was the first to calculate relaxation across different valley donor states. This 
body of work was the basis of some of the first proposals for electron spin qubits in silicon as a basis for a quantum 
computer [15J. This reinvigoration of the field has led to the extension of these original relaxation theories to new 
nanostructures like donors and quantum dots in strained silicon [3l[T6j[T7] and in III-V quantum dots [T8l H9j. 

The region of qubit interest in our case refers to a spin qubit with finite magnetic field, well below any degeneracy 
with higher orbital or valley states. We can summarize the key results of this paper and prior work on qubit relaxation 
relevant to qubit and quantum computer design in silicon quantum dots as follows: 

1. We are concerned with the lifetimes of excited states of a O-dimensional (OD) localized electron in silicon. The 
relaxation of ID and 2D mobile electron spins has been investigated for spintronics applications. This has led 
to some misconceptions about electron spin relaxation in ID and 2D vs. OD. They are in fact very different. For 
mobile electrons, scattering plays the key role [20) |2T], and spin rotation between or during scattering events 
is the driver of loss of spin memory (D'yakonov-Perel and Elliott-Yafet effects [22] ). This normally leads to 
spin lifetimes on the order of microseconds in silicon [2T| [23]. Since scattering is not an issue in OD, this is an 
unjustified worry for spin qubits for silicon quantum computers. However, spin lifetime measurements in silicon 
quantum wells can help determine relevant spin-orbit coupling parameters needed for quantum spin relaxation 
calculations. 

2. Spontaneous emission of a phonon is the dominant mechanism determining the spin- flip time, T\ (7\ being the 
characteristic time for true energy relaxation to the environment via a phonon or photon), at temperatures 
around 100 mK. Photon emission is negligible because of the much lower density of final photon states. Phonon 
emission accompanied by a spin-flip occurs due to spin-orbit mixing within the crystal. Of the "bulk" Roth 
[9] and Hasagawa [11] mechanisms that are relevant for donors in unstrained silicon, the "valley repopulation" 
mechanism disappears with increasing [001] strain as is common in silicon quantum wells [16J. The "one- valley" 
mechanism goes to zero if the magnetic field is parallel to one of the three crystallographic axes, and goes as 
B 5 for other directions \17\ . 

3. The effect of germanium in a SiGe QD heterostructure is not significantly detrimental to relaxation times for 
growth-typical Ge concentrations [16J (in the virtual crystal approximation). 

4. The spin-orbit coupling (SOC) in lateral Si quantum dots comes predominantly from structural inversion asym- 
metry (SIA) and symmetry-breaking due to interface effects leading to both Rashba-like and Dresselhaus-like 
SOC terms [22} [24 -26] . (The relative magnitude of Dresselhaus-like and Rashba-like SOC in silicon quantum 
wells or dots has yet to be verified experimentally, let alone systematically across samples.) Each term gives a 
characteristic magnetic field anisotropy in T\. Overall, 1/Ti is proportional to the seventh power of the magnetic 
field B |3J. This contrasts sharply with GaAs quantum dots and Si donor states. For these two cases 1/Ti oc B b 
(though for different reasons) . The ratio of the Rashba-like and Dresselhaus-like terms are expected to be sample 
dependent since in silicon they are due solely to interface effects (silicon, unlike III-Vs, has no bulk-inversion 
asymmetry (BIA)). 

5. Direct coherent rotations due to nearby spins in the semiconductor are possible. 

(a) At zero and low magnetic fields, direct dipole-dipole magnetic coupling with the central electron qubit and 
other electrons in the environment can occur. These rotations are technically coherent processes, but they 
result in spin flips that change the longitudinal component of the qubit magnetization and thus appear like 
T\ processes. These processes do not depend on B and can set upper limits on observed T\ times. The 
strength of this interaction is reduced as the inhomogeneity of the electron line widths increases, though 
even one electron spin exactly at resonance 200 nm away can cause 200 ms effective T\ lifetimes. A full 
theory is beyond the scope of this paper. 

(b) At zero fields, direct electron qubit - nuclei flip-flops are possible leading to Ti-like rotation. At finite fields 
this mechanism is suppressed due to the mismatch of the electrons and nuclei respective g- factors. 

(c) These mechanisms may in some cases be corrected via spin-echo techniques or suppressed by freezing out 
the background spins 

6. Other mechanisms for longitudinal spin relaxation, Ti, such as hyperfine coupling to Si 29 nuclei in natural Si, 
1/f noise, and Johnson noise, are estimated to be small in the parameter ranges considered here, though further 
work is required to verify this. The magnetic field dependence of these effects allows them to be distinguished 
from spin-phonon coupling. 
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7. The phonon bottleneck effect operates rather weakly in the parameter regime of interest for quantum dot 
applications. It can be calculated in a theory that goes beyond the electric dipole approximation; the result is 
only a slight enhancement of T\. 

8. Orbital relaxation in strained silicon is much faster than in bulk silicon in some important cases. Surprisingly, 
the rate of spontaneous decay from the first excited orbital state in silicon quantum dots can be comparable 
to that of GaAs quantum dots, which are commonly expected to relax more efficiently due to that crystal's 
piezoelectric nature. This has important implications for optically-induced motional spin-charge transduction 
(readout) and optical pumping (initialization) of spin qubits [27j, making the former harder and the latter easier, 
as well as for excited state spectroscopy. 

9. Excited valley state relaxation can be slow in silicon due to a small matrix element connecting valley states of 
different symmetry; this leads to the hope of valley qubits [4, 28 , 29J. The phonon emission rate has a maximum 
at the Umklapp phonon energy (11 meV and 23 meV for transverse and longitudinal phonons in silicon) that 
connects valley minima from one Brillouin zone to its nearest neighbor [13J. This leads to — at the longest — 
nanosecond relaxation times [30J in P donors with their large valley splitting (~10 meV), but in quantum dots 
is suppressed due to a large energy mismatch. However, the valley index in general cannot be considered a 
good quantum number in quantum dots due to large valley-orbit mixing [31]. For the same reason, "valley 
relaxation" in realistic devices can be dominated by orbital relaxation due to mixing with nearby orbital levels. 
The situation is different and more favorable in some cases such as for Li donors [H 130]. 

10. Non-ideal interfaces in silicon quantum dots may effect the relaxation processes; these effects are microscopic in 
origin [3TU35] and are not considered quantitatively in this paper. 

This content is arranged as follows. We begin with an introduction to the single electron states in silicon quantum 
dots typical of heterostructures used for qubits today. We follow with a discussion of the electron-phonon interaction, 
the dominant relaxation mechanism in these devices. We then use that theory to calculate the orbital relaxation (no 
spin flip) of low-lying excited states to states having the same valley index. We do this first using the electric-dipole 
(ED) approximation and then including all multipoles. This gives useful numbers for excited orbital state lifetimes 
as well as a quantitative idea of the extent of validity of the ED approximation in further calculations. Then the 
main subject of this paper is tackled, namely the spin flip mechanisms relevant to quantum dots. In this section, we 
review and adapt the previously known "bulk" spin-flip mechanisms to the quantum dot case. We then consider new 
mechanisms due to structural inversion asymmetry which give the dominant spin relaxation mechanisms for most of 
the magnetic field range. This is followed by comparisons with spin relaxation mechanisms (noise, nuclei) and valley 
relaxation, which completes the narrative of excited lifetimes in single electron quantum dots relevant to quantum 
computing. Where possible we compare the results for quantum dots with those for P donors in silicon and GaAs 
quantum dots, both of which are relevant reference systems. 

The final section summarizes the relation of theory to experiment. It is difficult to make sharp predictions for 
the absolute magnitude of T\ because of strong dependences on quantities that are usually somewhat uncertain, 
particularly the dot size, as well as a reliance on bulk material constants which may vary in nanostructures. We 
show how to overcome this problem by combining measurements of qubit spin lifetimes with measurements of excited 
orbital state energies and lifetimes. It is precisely for this reason that we deal in such detail with the excited states. 

II. SILICON QUANTUM DOT STATES 

This work is concerned mainly with lateral quantum dots formed by heterostructure confinement in the growth 
(z) direction and lateral confinement by metallic top gates. Figure 1 shows some of the heterostructure choices 
possible in constructing these dots, from modulation doped two-dimensional electron gas (2DEG) structures depleted 
by top gates to accumulation mode inversion layer devices in MOSFET-like structures. We will concern ourselves 
here predominately with the SiGe QW QD case, although our considerations should carry over to Si MOSFET dot 
structures as well. 

In semiconductors, the electron wave functions are superpositions of Bloch states at the bottom of the conduction 
band (CB), so the indirect band-gap, "many- valley" nature of silicon (as opposed to a single T- valley in GaAs) takes 
on great importance. In a biaxially-strained silicon QW, the number of states is doubled relative to GaAs, but reduced 
from the 6-fold degeneracy of electrons in the bulk. The conduction band (CB) minima located at k = (0,0, ±fco) 
with fco = O.S5k max have band energies lower than the minima at k = (±fco, 0, 0) and k = (0, ±fco? 0) by about 0.1 to 
0.15 eV at typical strain values (20-30% Ge in the barriers) [34J. Thus the k = (0,0, ±fco) valleys are populated but 
not the k = (0, ±fco, 0) an d k = (±fco? 0? 0) valleys [35]. A similar splitting is at work in Si MOSFET- type structures, 
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though here the physical origin of the lifting of the degeneracy is due to anisotropy of the silicon effective mass (and 
in some cases local strain). In the absence of magnetic fields and valley-splitting effects, the electronic ground state 
in this system is four-fold degenerate (spin and valley). When potentials or boundary conditions that mix the two 
valleys are present, as they always will be to some extent, there will be excited valley states corresponding to different 
linear combinations of the valley minima. So each valley state has its own identical set of orbital and spin states and 
an additional quantum number is needed to specify which valley state the electron occupies. 

A magnetic field splits the degeneracy of the spin states. The valley degeneracy is split by the hard confinement of 
the QW interfaces (or the impurity potential in that case) and is influenced by a number of parameters, especially the 
magnitude of the electric field in the growth direction and the sharpness of the confining potential. For the moment, 
we assume for simplicity that any static magnetic field is small or directed in the plane of the QW (so that the orbital 
functions are unperturbed) and that the well walls, located at z = {0,d}, are smooth. In these circumstances, we 
may write 

^=F£(x,y)F%\z) [a«u +z (r)e ife ^ + a« W -,(r)e-* " 

Fm (z) is the envelope function obtained by solving the confinement problem in the effective mass approximation for 
the m-th orbital, a±) z are the coefficients weighting the two valleys for the i-th valley state (i = 1, 2 for strained silicon, 
i = 1,2,3,4,5,6 for bulk silicon), and u± z (r) are the lattice-periodic parts of the Bloch functions, u 3 ;(r) exp(zkj • r), 

at the conduction band minimum kj. a^~^ = oS^J = l/y/2 for the symmetric valley ("sin-like") state and oi£~~^ = 

— o^lJ = 1/a/2 for the antisymmetric valley ("cosine-like") state. It is often convenient to expand the Bloch functions 
into a sum, 

u(r) = ^C G exp[ir • G], 

G 

where Cg weight the Fourier components of expansion (independent of r) and G are the reciprocal lattice vectors. 
The wave functions for a more realistic device, calculated in the tight-binding theory of Ref.[36j, are shown in Figure 
[8] (where the "Kohn-Luttinger" oscillations are evident but the lattice periodic oscillations are not included). A donor 
vs. dot energy level comparison is shown in Figure [2] 

Until now we have only concerned ourselves with "perfect" (completely flat) interfaces. In these cases the valley and 
orbital states are well defined - valley index being a generally good quantum number - much like the isolated donor 
case. In reality most interfaces are imperfect; they have alloy disorder (Si vs. Ge atoms), steps due to miscut, steps 
due to growth layer formation, and even interface states and traps (especially with respect to Si/Si02 interfaces). 
This leads to mixing of the valley and orbital wave functions as well as diminishing of the valley splitting (due to 
interference) [3Tjl33]. and as such, wave functions and splittings that vary from device to device and dot to dot. These 
microscopic variations are not considered in this paper. Fortunately, these effects are often not large corrections to 
the calculations below, as the experimentally accessible energy splittings come into the equations at a much higher 
power than the relevant (experimentally inaccessible) matrix elements. As we go we will point out differences from 
our theory for the ideal and likely situations, focusing on experimentally accessible signatures. 



III. ELECTRON-PHONON INTERACTION 



Phonons are quanta of lattice vibrations, that is, mechanical motion that gives time-dependent stress. This alters 
the band structure by shifting band energies and lifting degeneracies [35j [37] . It is typically assumed that this effect 
does not change the band curvature (effective masses) but does shift the energy states of interest. The shift in energy 
of the band edge per unit elastic strain is called the deformation potential and is common to all semiconductors and 
solids. In polar crystals, distortion of the lattice can also create large internal electric fields which affect the electron. 
This is the piezoelectric interaction. Ionic crystals like GaAs suffer from piezo-phonons, which are often very efficient 
at electron scattering; silicon, being non-polar and centrosymmetric, has none. While optical strain also exists in 
materials with two atoms per unit cell, such as silicon, optical phonons in silicon have a narrow bandwidth centered at 
a much higher energy than a quantum computer would operate. We briefly review the theory of the electron-phonon 
deformation potential interaction outlined by Herring and Vogt [38] . 

The energy shift of a non-degenerate band edge due to strain is given by 

HeL = J2u a0 E%, (1) 

a,/3 
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Figure 1: Different types of quantum dot heterostructures relevant to the theory in this paper. 



Constant 


Value 


e 


1.6 x 10" 19 C 


h 


1.05 x 10" 34 J s 


c 


3 x 10 8 m/s 


eo 


8.85 x 10" 12 C 2 /Nm 2 


£Si 


11.8 


All (Si) 


1.999 


SJ-(Si) 


1.998 


H„(Si) 


8.77 eV=1.4xl0" 18 J 


H d (Si) 


5 eV=8xl(T 19 J 


p(Si) 


2330 kg/m 3 


wj(Si) 


9330 m/s 


wt(Si) 


5420 m/s 


ao(Si) 


0.543 nm 


kmax (Si) 


27r/ao 




1.38xl0" 23 J/K 



Table I: Physical constants and materials parameters for bulk silicon. 



where 



a/3 



dui 
dxi 



duj 
dxi 



is the strain tensor and SW is the deformation potential tensor for the ith silicon CB valley. Since our electron 
is confined to massively strained [001] silicon, we need only include the (0,0, ±fco) valleys which have deformation 
potential tensors given by 
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sin# 



Table II: Polarization components. 



where relates to pure dilatation and S n is associated with shear strains. is a unit vector in the direction of 
the ith valley. For silicon under compressive stress along [001], opposing valleys move in energy identically. In order 
to associate phonon modes with strain, we can expand the unit cell displacement u(r) in plane waves, 

u(r) = J2 [ e fa> A ) + e* (q, A) a qX e~^] . 

qX 

a q \ destroys a phonon with wavevector q and polarization A (2 transverse and 1 longitudinal; see details in Table 2) 
of a phonon; e (q, A) is its unit displacement vector. This results in a strain tensor due to a phonon of 



U (q, A) a/3 = - [(e a (q, A) qp + (q, A) <? a )a* A e * q r + (e a (q, A) qp + (q, A) q a )a q \e ZCL ' r ] . (3) 
The operators a q \ and a* A have matrix elements 



(n qA - 1| a qX \n qX ) = Jhn qX /2M c u qX , 



(n qX + 1| a* A |n qA ) = y (^ q A + 1) /2M c cj qA , 

M c is the mass of the crystal and n q \ = 1/ {e huJc ^ x ^ kT — l) is the phonon occupation number of the mode with wave 
number q and polarization A. The complete electron-phonon Hamiltonian must be summed over phonon modes and 
polarizations. For a [001] strained-silicon quantum well, it can be written succinctly as 



H zp = iq t a <?A e * q ' r + a q\ e%Ci ' r ] x i-d^x (q, A) q x + E d e y (q, A) q y + (S d + E u ) e z (q, A) q z ) . 



(4) 



A=l q 



We are especially concerned with the anisotropic effects due to the massive strain of the system in question. As can 
be seen from Eq. [2] and the deformation constant values in Table |TJ the shift in energy of a specific valley due to an 
acoustic phonon is very anisotropic. In the case of bulk Si, the six conduction band minima are equidistant from the 
T-point and thus form an isotropic response to phonon deformations. This means essentially that transverse phonons 
will not contribute to the relaxation times for intervalley transitions of the same symmetry (that is, otfi) = ol$) for 
initial state i and final state j). Another way to see this is to consider the electron-lattice matrix element between 
different plane wave states at the same minimum (Equation 3.29 of Ref. [llj), 



H, 



eL 



(0 



qt 



ie(q,t).3W-q 



J mn 



(q) 



C.C. 



where /mn(q) = / Frn (r)e ZCL ' r Fn^ dr. We have used the polarization index t to indicate a transverse phonon. The 
matrix element between two dot wave functions is then 



m t m 



H, 



eL 



qt 



ie (q,t). VaWaW 



3W • q 



/^(q) + c.c. 



It's easy to see from the above equation that if proportional to the identity matrix (assuming 

<^m = <^n = 1), then the transverse phonon matrix elements must be zero since e t i(q,tl) _L e^(q, t2) _L q. The point 
is that in strained silicon, only the ±z minima are occupied so unlike the bulk silicon case, transverse phonons will 
contribute. This turns out to be very important in relaxation calculations, as will be seen below. 
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Figure 2: Energy level diagram for P donors in Si (left) and Si quantum dots (right). Note the reduction in energy scale from 
left to right. In both cases, the electron exists as a superposition of different conduction valley minima within the crystal. 
Different low-lying dot levels (all s-like) result from the sharp donor potential and different amplitudes in the six valleys shown. 
Low-lying dot states arise from different orbital wave functions with a further valley splitting (far right) between symmetric and 
anti-symmetric combinations of the two valleys, zLz. The circled numbers represent relevant relaxation processes: (1) orbital 
relaxation across the first energy gap; (2) spin relaxation of the electron qubit; and (3) valley relaxation. Note that the valley 
splitting can vary from to meV and may fall within the energy levels of the low- lying Zeeman and orbital states. 



IV. ORBITAL RELAXATION IN STRAINED SI QUANTUM DOTS 



The line widths and characteristic behavior (dependence on magnetic field, etc.) of the lowest lying excited states 
in quantum dot systems are usually very relevant for characterizing a spin-based qubit. We begin by considering the 
relaxation of an excited state that involves no spin flip and takes place within the same valley state (assuming for 
now no valley-orbit mixing) . Spontaneous emission of a single acoustic phonon is the dominant relaxation mechanism 
for excited electronic states in Si at QC temperatures (<100 mK) (an ansatz guided by empirical evidence for silicon 
donors [6[ [39] and GaAs quantum dots |40j). The deformation potential approach [38J is easily modified to include 
strain [16J. The phonon-induced energy shift is very anisotropic for a silicon conduction band valley. Because of this 
the results for strained silicon are quite different from those of bulk silicon [8 J . 

In the bulk, the six conduction band minima are equidistant from the T-point and thus form an isotropic response 
to phonon deformations (specifically the case where ctj = 1 for all i across both states). So the angular integral over 
transverse phonons averages to zero. For example, the transition from the ground, ls-like symmetric state to the 
2p-like, symmetric state in bulk silicon has no transverse phonon contribution |8J. The same transition in strained 
silicon does have such a contribution, because the cubic symmetry of the six minima has been broken by strain, and 
indeed it is the largest term. This greatly increases the relaxation rate since the inverse sound velocity comes into the 
rate equations with a very high power as we will now show. 

Fermi's Golden Rule, between states of arbitrary spin, 

\{nsv\H ep \ms'v')\ 2 S (E ph - E nS)Tns /) , (5) 

is the basis for our phonon relaxation rate calculations. Here, \nsv) is the state ipns of the electron in the dot on 
level n with spin state s and on valley manifold v, including the effect of a magnetic field. We assume an isotropic 
phonon spectrum such that the energy of the phonon is E p h = ^j q A, where cj q A = v\\q\ and ^a is the velocity 
of the mode A. Setting s = s' for orbital relaxation without a spin- flip and v = v' , E nSjfns f = E mn is the energy 
splitting between states m and n, H ep is the electron-phonon interaction of Eq. H] Summing over phonon modes, 
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using q = (sin cos (/>, sin sin 0, cos 0) and e/ = g _L e t i _L e t 2 (see Table [TT|) for the wave vector and polarization 
vectors and using the electric dipole (ED) approximation, e* q r ~ 1 + iq • r, we find for the phonon-induced relaxation 
rate, 



\E n 



h 6 7rp S i 



Mi r > 



y 



- xy 



Mi r > 



T z \ (n B (E mn ) + 1) : 



where 



T — 

xy 



35S2 + 14S d S tt + 3SS , 2~i 



210vJ 



105v 7 t ' 

355^ + 425,5, + 155^ g 
+ 35vJ' 



(6) 



(7) 



210vJ 



the matrix elements are M( mn ) = (F m |r|F n ), is the mass density of Si, E mn is the energy gap between orbital states, 
V£ and v t are the longitudinal and transverse sound velocities. The single electron envelope functions can be calculated 
by solving the Poisson and Schroedinger equations directly as in Ref. [lj or, as is normally done, by approximating the 
potential as a parabola, giving harmonic oscillator states defined by the fundamental energies E\$ = ftw Xjy . With the 



envelope functions F = (2 /it) 
is given by 



1 ^ 4 Xq 1 ^ 2 exp (— x 2 /x 2 ) ) and F\ = (2/7r) i/4 (2/ y /r x^)xexp (— x 2 /xq), the matrix element 



a/4 



Mi 10 ) 



m t E 



10 



(«) 



where 2^(x 2 ) = max{xo, yo} = V2h/\^m t A is the lateral size of the dot, L, and xo and yo are the dot sizes in the x 
and y directions. We define A = E\ — Eq. Thus, the orbital relaxation rate from the lowest orbital state for a slightly 
asymmetric (non-degenerate excited state), parabolic dot is 



25 2 



lObvJ h 2 7rpsim t 



(n B (A) + 1) 



(9) 



(we have used the fact that v\ ~ 2v t to eliminate terms due to longitudinal phonons). Because ^ mn appears in the 
fifth power (general case) or fourth power (parabolic dots) in Equations [6] and [9] for orbital relaxation, an accurate 
value for E mn is much more important than equivalent accuracy in the wave function matrix elements. The dipole 
approximation is valid until roughly qL » 1 {L is the maximal linear size of the dot) when the relaxation rate starts 
to decrease due to phonon bottleneck effects. We treat this effect explicitly in the next section. 



V. PHONON BOTTLENECK EFFECT 



When the dot becomes very small [L comparable to a few interatomic spacings), the relaxation rate is reduced. 
This is due to the impossibility of satisfying simultaneously energy and momentum conservation during an electron- 
acoustic-phonon scattering event |41| . Mathematically, the phonon bottleneck effect is due to the fast oscillating 
exponential factor: emission of a phonon with wave vector q is unlikely when q > 2tt/L [42J. Table ?? charts this 
transition for parabolic Si quantum dots. Note that in a lateral quantum dot, unlike excitonic quantum dots, Auger 
processes and electron-hole scattering do not play a role in negating phonon bottleneck. Assuming parabolic dots, 
only dots with fundamental energies A « 1 — 2 meV are small enough (and virtually impossible to construct with 
laterally-gated devices) for phonon-bottleneck effects to make a significant impact on increasing the orbital relaxation 
times. 

We outline the calculation for Si, including the valley effects, in Appendix |IX A| The chief difficulty is to include all 
multipole moments. The result is that the orbital relaxation rate for a parabolic dot (in all three dimensions) from 
its first orbital, excited state to its ground state is given by 



Tt act ={n B (A) 



exp 
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A (meV) 


L (nm) 


0.05 


127 


0.1 


90 


0.2 


63 


0.3 


52 


0.4 


45 


0.5 


40 


1 


28 


2 


20 


3 


16 


8 


10 


10 


9 



A (nm) - Si(v t ) 1/F ED (Eq. 6) l/Y exact (Eq. 10) 



447 
223 
112 
75 
56 
47 
22 
11 
7.4 
2.8 
2.2 



2.5 x 10" 7 s 

1.6 x 10" 8 s 

9.7 x 10~ 10 s 
1.9 x 10" 10 s 
6.0 x 10" 11 s 

2.5 x 10" 11 s 

1.6 x 10" 12 s 

9.7 x 10~ 14 s 
1.9 x 10" 14 s 

3.8 x 10" 16 s 
1.6 x 10" 16 s 



4.0 x 10" 7 s 
3.8 x 10" 8 s 

5 x 10" 9 s 
1.8 x 10" 9 s 

9.5 x 10" 10 s 

6 x 10" 10 s 
1.4 x 10" 10 s 

3.8 x 10 -11 s 

1.9 x 10" 11 s 

9.6 x 10" 12 s 
1.4 x 10" 11 s 



Table III: Characteristic numbers for the phonon bottleneck effect in Si and orbital relaxation rates. A is the first excitation 
energy, L is the lateral dimension, and A is the wavelength of the transverse phonon with the resonant energy. 



where zo is the size of wave function in the growth direction and q&i = A/hvi, qAt = A/hvt. The A coefficients are 
defined by 



4 n) (9A.) = f i 1 x n exp (±{x% - zl)ql 8 x 2 



The results are plotted in Figure [3j Note that the exact expression, Eq. [TUJ reduces to the dipole approximate 
expression, Eq. pi when q\ = as expected. For large A, the exact solution for the orbital relaxation rate Ti2 begins 
to diverge from the electric-dipole approximation early on and never falls below a picosecond or so. Despite this, the 
electric-dipole approximation holds well for small energy gaps, 0.1 — 2 meV, where a quantum computer in silicon will 
most likely operate. This figure shows that there is no significant benefit in going beyond a few meV. Only around 10 
meV does the relaxation time start to increase, but this a relatively small effect. It does demonstrate, however, that 
the phonon bottleneck effect may be experimentally observable in these systems and, more importantly, that our use 
of the electric-dipole approximation gives results below for the spin-flip times that may be considered a lower bound 
on the maximum possible time T\. 

Table [TTT] compares the output of Eq. [6] and Eq. |10| for pure orbital relaxation in lateral silicon quantum dots. 



One can compare these results to those for GaAs quantum dots [43], where typical values would be 1/r = 10 -8 s for 
A = 1 meV and piezo-phonons dominate. Excited orbital states in strained silicon typically relax in nanoseconds or 
faster, corresponding to a level broadening of a micro-eV or wider. It is possible that a low- lying excited valley state 
(of the same spin direction) may be closer in energy than the orbital level. For bulk silicon, theory and experiment 
have found characteristic relaxation times for the 2p-ls transition of ~200 ps [44J. 



VI. SPIN-FLIP (Ti) TIMES 

Our expressions for orbital relaxation in strained silicon can be extended to the case of a spin-flip transition due 
to spin-orbit coupling (SOC) in a QW. We expect that relaxation via a phonon is the dominant cooling mechanism, 
in this case mediated by SOC which mixes pure spin states. This is known to be the case for donor-bound spins in 
bulk [6J. Structural inversion asymmetry has traditionally been thought to be the dominant source of SOC in silicon 
quantum wells due to the large electric field common to modulation-doped or top-gated SiGe heterostructures. The 
nature of this SOC has been well described elsewhere [21] I24H26) [45] : it leads to a Rashba term in the Hamiltonian of 
the electrons. However, interface effects that break the inversion asymmetry can also lead to a generalized Dresselhaus- 
like term [25]. Surprisingly, this can lead to effects of similar or even greater magnitude than the Rashba term [24U26] . 
We discuss this further below. Here we only note that the zero-field energy level splittings caused by SOC are small, 
of the order of /ieV, which validates our use of a perturbation theory that uses zero-order electron wave functions and 
energy levels taken from the SOC-free Hamiltonians. 

Until the appearance of Refs. OUT], there was no finite spin-flip time prediction for lateral silicon quantum dots 
when the external field B is parallel to z. Previous theories for T\ in silicon had been based on the two dominant 
mechanisms relevant to P:Si donors: the "valley-repopulation" mechanism (bulk SOC mixing with the six nearby 
ls-like states) and the "one- valley mechanism" (bulk SOC mixing with continuum states) [9HTT]. Both mechanisms 
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i/r (s) 1/r (s) 




Figure 3: Relaxation rates, 1/r, for the transition from the first excited orbital state to the ground state by emission of a 
phonon versus orbital energy splitting, A. The exact results including all multipole contributions are given by the blue solid 
lines while the electric dipole approximation results are given by the purple dashed lines. 

are independent of the size and shape of the localized electron wave function. We showed rigorously in Ref. [16] that 
the former becomes negligible with [001] strain. The latter is slightly modified with strain and goes to zero for certain 
directions of the static magnetic field, particularly the [001] direction (the most relevant to QC), for both one and 
two-phonon processes [17] . We review these bulk mechanisms here as they are relevant for donor qubits and in some 
cases may be seen as residual T\ mechanisms at low magnetic fields in dots. Then we will derive the spin relaxation 
times for dots in strained structures due to inversion asymmetry-based SOC leading to a T\ that is finite for B parallel 
to z. 



A. Bulk spin-flip mechanisms 

For donor states in bulk silicon, Roth and Hasegawa [9HTT] identified two mechanisms that have been confirmed 
experimentally up to T =2 K in P donor spins [12J. The spin relaxation in the Roth-Hasegawa picture is due to a 
modulation of the system's g-factor by acoustic phonons. Both mechanisms are direct single-phonon processes. The 
g-tensors for a given conduction band state can be written as a sum over the g-tensors at each conduction band 
minimum, 

g = ^aigi, 

where |c^| 2 is the squared amplitude ("population" in the early literature) of the single electron wave function at the 
zth valley and g w = g±S a(3 + (g\\ - g±)K$K%\ 
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There are two mechanisms leading to spin flip. The "valley-repopulation" mechanism is due to mixing between 
the symmetric ground state of the donor electron and the split-off doublet state where a phonon changes the a[s. 
The "one-valley" mechanism is due to phonon-induced modulation of the themselves and subsequent mixing with 
nearby conduction bands which are coupled through an inter-band deformation potential. The two mechanisms are 
of the same order of magnitude in the bulk case and complementarily explain the angular magnetic field dependence 
of T\ in those systems [T2j . 

Ref. [16J showed how the valley-repopulation contribution to the spin-flip becomes negligible with increasing [001] 
compressive strain as is inherent in a silicon quantum dot. This can be seen easily qualitatively. Consider first a 
potential with spherical symmetry. The population amplitudes describing the lowest six conduction states are given 
by [12: 

Singlet : an = -^(1, 1, 1, 1, 1, 1) 
Doublet: a 2i = ^(-1, -1, -1, -1, 2, 2) 

a 22 = ^(l, 1,-1,-1,0, 0) 
Triplet: a 3 i = ^(1, -1, 0, 0, 0, 0) 

^32 = ^(0, 0, 1,-1,0, 0) 

^33 = ^(0, 0, 0, 0, 1,-1), 

in the valley basis (x,—x, y,—y, z,—z). The six states are split even at zero strain by non-spherical central-cell 
corrections in the donor case (the sharp potential of the donor) and by the interfaces in the quantum dot case. SOC 
represented by the anisotropic g-factor mixes the singlet ground state only with one of the doublet states [16j. In 
QDs the strong compressive strain in the z direction causes a large relative splitting between the six valley states; the 
result is that only the an and as3 states will be populated. These symmetric and antisymmetric valley states are not 
mixed by the SOC which results in a vanishing matrix element. So, in the quantum dot limit (±z valleys populated), 
the valley-repopulation contribution to the spin-flip rate becomes negligible. Note that for this mechanism we only 
consider mixing to the six lowest states of the donor, all of which have orbital s-like character. The 2p states in a 
donor are typically 30 meV (bulk) to 3 meV ([001] strain) away. Mixing with these states will be considered separately 
below. 

The one-valley mechanism, however, is relevant to QDs. Roth showed that in bulk silicon, the contribution from 
mixing with nearby bands is described by a Hamiltonian 

H b t k - vaUey = Ap(U xy (a x H y + a y H x ) + c.p.), 

where c.p. stands for cyclic permutations. Group theoretical considerations and perturbation theory lead to the 
conclusion that the dominant contribution to A comes from mixing with the nearby A 2 and A5 bands and is given by 



A = 



2i/3 (A 2/ |p,|A 2/ )(A 2/ |^|A 1 ) 



3m E\ 2 ,E Xh 
x {(Axl^lA* ><A*|/gA 2 ,> + (A 1 |^|AD(Alb £C |A 2 ,)} 



xA^|^(A 2 ,|^|A 1 ), 
^12 



where h = VV x P is the usual crystal spin-orbit vector, are the energy gaps to the relevant bands, and D is the 
inter-band deformation potential. H^ l ^_ valley represents a sum over the six minima in the bulk case, but in a lateral 
quantum dot the dominant contribution comes only from the ±2 minima and was determined by Glavin and Kim [17] 
to be 

H^ne-valley = Af3U xy (<J x H y + <J y H x ). 

The constant A was experimentally determined by Wilson and Feher [12J as A = 0.44. The T\ time due to H®^_ valley 
can be readily calculated and is given by [17] 

1 2v i ^ A/* Pt ( B \ ^ 
ji— = ~^~~o — r f tt^ I (1 + 2n B (^/i5))sin 2 (9(cos 2 20 + cos 2 (9sin 2 20), 

T one-valley $ g 2 py 5 ^ 2 ixh J V V Y ^' 

where (#,(/>) define the angle of the magnetic field relative to [001]. n# is the Bose function. It is evident that this 
equation produces infinite relaxation times if the magnetic field points along the [001] or [011] axes. Figure [4] plots 
the one-valley relaxation rate as a function of magnetic field direction. We particularly stress that the B 5 dependence 
of 1/Ti is the characteristic signature of this mechanism and that this rate is very small in general. 
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Figure 4: Dependence of the spin relaxation rate on the angle of the external magnetic field for the one- valley mechanism in 
strained silicon quantum dots (or donors in strained silicon). 



Now we will consider SOC that comes from the dot structure itself. The spin-orbit Hamiltonian for a two dimensional 
electron is conventionally written as 



where k x , k y are the in-plane wave vector components; a and (3 are the strengths of the so-called Rashba and Dressel- 
haus spin-orbit terms. This, combined with electron-phonon coupling, can also produce spin relaxation, and the effect 
(with a = 0) has been computed for GaAs QDs [18J. In Si QDs, there is SIA SOC that comes from the fact that the 
mirror symmetry z — z is broken. This occurs in SiGe/Si/SiGe heterostructures and in MOSFET-type QDs either 
by modulation doping or by a top-gate induced electric field [2T] , In contrast, Dresselhaus SOC has traditionally been 
assumed to be absent in Si structures because bulk Si has inversion symmetry. It was recently shown [25 , 26J that this 
is not the case. The breaking of inversion symmetry by the interfaces gives a non-zero Dresselhaus-like term that can 
be surprisingly large [26J. This has yet to be verified experimentally. Often, experimental measures of spin relaxation 
involve terms proportional to a 2 + (3 2 so the terms are hard to verify independently. Therefore we will keep both 
terms in Hso and compute the spin relaxation that comes from these asymmetry-induced effects. The results differ 
from those of the GaAs quantum dot analog due to the many- valley nature of silicon and the dominance of acoustic 
over piezo-phonons in silicon. 

The orbital energy level splittings are much reduced in quantum dots relative to donors (see Figure 1) because of 
the more shallow potential. It is thus relatively easy to make the Zeeman splitting larger than the orbital splitting. 
However, for quantum computing the likely situation is for the Zeeman splitting to be less than the orbital splitting 
to maintain a good qubit manifold. Here we will consider only this case where the magnetic field splitting is small 
compared to the orbital splitting. We comment on this approximation further below. In Si, the Hamiltonian for the 
electron-phonon matrix element is, for v — v'\ 



where \ns) denotes a state with an electron in the n th level of the dot with spin s. Again, we consider relaxation 
processes within the same valley state. A phonon with wave vector q and polarization A is absorbed or emitted 
depending on whether E n > Em or E n < E m . s =t, I is the spin projection on the z-axis, defined to be along the 



B. Dot-specific mechanisms 



(ms \H ep \ ns') qA = i$s,s> (m | [a* A e iq ' r + a qA e iq ' r ] | n) 

q [(E d e x (q, A) q x + S d e y (q, A) q y + (H d + S u ) e z (q, A) q z )\ 
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e-^/ 2 cos9/2 



e 4 */ 2 sine/2 / ' 



external applied field B = B (sin O cos <£, sin 6 sin <£, cos G) . Hence |t) 
The matrix elements of H so are 

(ms\H so \ns f ) = (ms\[a (cr x ky - cr y k x ) + f3 (cr x k x - cr y k y )]\ ns') 



etc. 



imtEry 



where x mn = (m\x\n) is the dipole matrix element for the dot states, <5\j^ = — cosy? cos # — i sirup and <r^ = 
— sin^cos^ + icosip, and where we have used the trick (p x )kn = imEk n %kn/fy we use units with ft = 1. H so causes 
the eigenstates to be mixtures of up and down spin states. For example, if the unperturbed orbital ground states 
|0 t) (0) and |0 |) (0) are perturbed by H so , the new eigenstates |0 t) (1) and |0 are 

|0 t) (1) « |0 t) (0) + im t J2 C 1 + 9^BB/E m ) ((ay mn + /3x mn ) <rf - (ax mn + /fy mn ) ) \m |) , 

ra/0 

|0 « |0 |) (0) + im £ ^ (1 - g/n B B/E m ) ((ay mn + /3x mn ) cr^ - (ax mn + /?2/ mn ) cr^) |m t) , 

ra/0 



^(l T ff + ...). The (1) 



where m t is the transverse mass and we have expanded around g\iB with g k ± g ^B 

superscript indicates first order in \H so \ /E m - the spin-orbit splitting compared to the orbital excitation energies. It 
is important to compute the gfisB / 1 E m correction for reasons that will soon become apparent. 
Our interest is in the matrix element 

<0T|ff ep |0|} qA , 

where the qX subscript indicates that there is a phonon in the final state. We find 

<0t|tf e p|0|) qA = ^ <(0)|a qA |(l))(z 

[(E d e x (A) q x + E d e y (A) q y + (E d + E u ) e z (A) q z )\ x 

(0 t| (0) - im t + gi^BB/E m ) ((ay 0m + fem) ^ - (ax 0m + Pyom) erf)* (m i\ 

|0 |) (0) + im £ (1 - gUBB/Em) ((ay 0m + ^ 0m ) - (ax 0m + /3y 0m ) cr^) \m t) 

ra/0 

and we make the electric dipole approximation e* q r w 1 + iq • r, which gives 

(0 t \H ep \ Dfr = -2im t ^ r m0 ? (O^Om + ^Om) ^ ~ (^0m + fem) ^) 



=x,y,z ra^O 



2M c cj g -* A 



[(S d e x (A) $x + S d e y (A) q y + (S d + S u ) e 2 (A) ^)] ^. 



Because (Hso)kL = — (^scOtl? the overall matrix element is reduced by roughly gfiB/E m . This is the manifestation 
of the so-called Van Vleck cancellation. 

We do the thermal average over phonon states and apply Fermi's Golden Rule. This yields 



±r = 2tt J2 [1 + (u^)] |(0 t \H ep \ i)^ S (gfi B B - v x q) 



A A m,n/0 

^ (<f • r m0 ) (g • r n0 ) [(S d e x (A) ^ + (A) ^ + (E d + E u ) e z (A) £,)] 2 , 
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where 

Sm = ((Otyom + PxQm) (7^ - (ax 0m + PVOm) ■ 

The integral is over the directions of q = q (sin cos 0, sin # sin 0, cos #) . (6 and are not the same as G and <I>, which 
give the directions of the magnetic field.) 

Following Ref. [19] we now define the dot polarization tensor 



= -2e 2 iXi)m ° ^ 



Ora 



where the sum is over all the orbital states. The opposite valley states are not included in the sum since the intervalley 
electron-phonon coupling is assumed to be small (this could be different in non-ideal interfaces). It is also reasonable 
to neglect z m o, since the spatial extent of the wave function in the growth direction is small compared to £o,2/o ( at 
least by a factor of 10). The result (with h restored), is 

^ = T xy -g {g^ B B) 7 + 2n B (g» B B)} x 

1 1 Till psi 6 

{ [(a + f3 ) (£xx£xx H~ ^yx^yx H~ £,xy£,xy H~ £,yy£,yy) H~ 2a/3 (^^^^ + ^yx^yy H~ £xy£xx H~ 

3 1 N 
- + - cos 26 

4 4 

"I - ) (^xs^a; H~ ^yx^yx ^xy^xy ^yy^yy) 

x sin 2 6 cos 2$ 

~^ 2 ^ ) (^> XX ^> X V ^>U x ^>yy ^>xy£>xx H~ ^yy^yx) 2tt/3 (£##£## H~ ^yx^yx H~ ^xy^xy H~ £2/2/ £2/2/)] 

x sin 2 6 sin 2$}. (11) 

Note that we include a factor of two in the phonon population multiplier, 1 + 2n&, so to satisfy the traditional 
definition of T\ where both relaxation and excitation are possible (although at very low temperature this term goes 
to one). The most striking qualitative feature of this expression is the B 7 dependence |3J; this can be considered as 
the characteristic feature of dot-specific SOC and contrasts with the B 5 dependence of bulk SOC, as well as GaAs 
quantum dots (which are dominated by piezo-phonon relaxation in energy regimes of interest [19J). In addition, there 
is field anisotropy. To understand this anisotropy note that the diagonal elements £ xx « £ yy are likely to dominate 
the off-diagonal elements ^ xy and £ yx . Examination of the expression then shows that the largest term in 1/Xi is 
proportional to (a 2 + /3 2 ) (3 + cos 20) /4. This does not vanish when B is along the z-axis (O = 0), again in contrast 
to the Roth-Hasegawa contributions. If one wishes to determine a and f3 individually, then the smaller contribution 
proportional to {a 2 — P 2 ) sin 2 O cos 2$ must be measured. It could be enhanced if £ xx is very different from £ yy which 
would be the case for a very elliptical dot. 

The polarization tensors (matrix elements) can be calculated numerically, or in the small magnetic field limit (where 
the first excitation energy is much less than the Landau energy), the zero B-field parabolic matrix elements can be 
used as a decent approximation. We can include the magnetic field in a circular dot explicitly if B\\z by utilizing the 
Fock-Darwin states [50]. In this case (00 1 x |01) = <00| y |10) = \ V2L B /2\ , where ficj is the fundamental energy of the 
dot, Lb = 2^^/ (m*Q), Q = \/oJq + cj^/4, uj c = eB/m* , £?oo,oi = fi^-, ^oo,io = and uj± = v^o + uj 2 /4:±uj c /2. 
Then, taking into account dipole selection rules (only transitions n± = n± ± 1 are allowed), 



£yy — e 



2 I L _B 

2 
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hjjj- 





aI1 d ^ X y — £y X — 0, 



In the 5 = limit, with L = Lb(B = 0), reduces to £ xx = % 2 and £ yy = % 2 for an elliptical dot. 



To estimate the overall quantitative magnitude we shall assume a circular dot with a parabolic potential (note also 
that for a circular dot i xx (B z ) = £ yy (B z ) = £ xx (0). Then with £ xx = £ yy = f and £ xy « £ yx = and 

1 _ 1 . 2 (« 2 + / ?2 )(^) 7 ( 3 + cos2e)[1 + 2ns( ^ )] (12) 
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Fig ure 5: Spin relaxation time, Ti, for an ideal, circular quantum dot as a fun ction of magnetic field as calculated from Eq. 
I12I for SOC constant value (which can vary greatly from device to device) of \J ol 2 + f3 2 = 4/xeV-nm (taken from a SiGe QW 
experiment [46 - see text) and 5||[110] (T\ for B\\z is a factor of 2 longer). Diagonal lines for spin relaxation times from bottom 
to top are for different orbital energy splittings (dot size gets smaller going up) of A = 0.02 meV, 0.05 meV, 0.1 meV, 0.2 meV, 
0.3 meV, 0.4 meV, 0.5 meV, 1 meV, 2 meV, and 8 meV, respectively. Theory is only appropriately compared to experiment 
below degeneracy points (given by intersection of vertical and diagonal lines). Any possible spin "hot spots" would occur at 
degeneracy between Zeeman and first excited state splitting (which is approximately 0.3-0.4 meV in these size dots). Points 
represent presently published experimental data for a SiGe quantum dot from HRL (squares) [37], a SiC>2 quantum dot (circles) 
from UCLA [48], and a SiGe quantum dot (triangle) from Wisconsin [49] . 



where we have used 



^ (13) 



i t 2 A 4 2 4 /L 4 



xy 



(for reference we find A = ) • The contribution of longitudinal (£) phonons is suppressed by roughly a fifth in T ; 
and is neglected. For reference, the spin relaxation rate can also be written in terms of the dipole matrix elements 
between Is and 2p, M^ 10 \ as (assuming mixing to one excited state): 



1 = 4 ^ 2 (a 2 +/3 2 ) (g^ B B) 7 \M 
T QD ~ 1Q5 *u nfppsivj A 2 



(3 + cos 26) [1 + 2n B (gfi B B)] . (14) 



The magnitudes of a and j3 are material system and device specific. Wilamowski et al. [46J have measured 
aw = 0.55/ a/2 x 10 -12 eV-cm = 4 fieV-nm via 2DEG spin relaxation in Si/SiGe quantum wells. On the other hand, 
for a SiGe/Si/SiGe well and a field of 10 7 V/m (roughly a factor of 2 larger E- field than typical SiGe QW QDs but 
about right for Si02 dots), Prada et al. [26J theoretically find that j3 > a and j3 = 5.77 /ieV-nm = 9.2 x 10 -34 J — m; 
note that they mention that the j3 term could decrease in (typical) heterostructure quantum wells with miscut, and 
that it will very from device to device. Calculations of the spin relaxation in 2DEGs [20j[2T] give \/t1 deg oc [a 2 + j3 2 ) 
assuming minimal cyclotron effects or B \\ (x,y); so the two results are consistent if we attribute the Wilamowski 
result as due to /3. Nestoklen et al. [51], also theoretically, find a value for /3 roughly 6 times smaller than Wilamowski 
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Figure 6: Schematic summary of the mechanisms behind the different field dependancies of the spin relaxation rate for GaAs 
dots, for which 1/Ti ~ B 5 , for Si donor states dots, for which 1/Ti ~ B 5 , and for ideal Si dots, for which 1/Ti ~ B 7 . 



et al. [46J (Pn — 15.6, a at = 5.2 /xeV-nm). Note that the measured line width in Ref. [46J does not depend on the 
in-plane orientation of the magnetic field, implying (at least for that device) that one SOC term dominates over the 
other |20J|21]. 

Figure [5] plots Equation 12 for B || [110] as a function of orbital energy splitting using the result of Wilamowski 
et al., aw = \J {ol 2 + P 2 )-, and substituting values from Table [tl Also shown are some recent experimental results 
[47H49] which follow the B 7 trend but generally show longer lifetimes. We defer our comparison to experiment to 
the Discussion section below. Note that our results are for the electric dipole approximation, shown in the previous 
section to be a good approximation below ~1 meV. Calculations that attempt to calculate the spin relaxation beyond 
the dipole approximation [52j [53] have shown the possibility of fast spin relaxation ("hot spot") when the orbital and 
Zeeman energies are degenerate: a result of orbital-spin level mixing. This mixing is contingent on the nature of the 
SOC mixing — Dresselhaus (no mixing) or Rashba (mixing). Unfortunately we do not know the nature of the SOC 
in these devices and indeed it may depend on microscopic details and very from device to device (even on the same 
chip) . An additional complication arises when the first excited state is a valley state or valley-like - a situation we will 
discuss later in the text. So a hot spot is not assured. Therefore, we mark these crossovers in energy in Figure ?? as a 
possible position of interesting physics (fast relaxation) which may tell us more about the nature of these states. Our 
results are relevant for the quantum computing situation where g\iB < A and likely a good approximation beyond 
the crossing point (with the cyclotron modified wave function incorporated). 



C. Electrical and Magnetic Noise 

Another possible mechanism that could limit T\ in Si/SiGe quantum dots is the electric and magnetic noise coming 
from trapped charges and other two-level systems, noise in the circuitry, thermal and quantum current fluctuations 
in nearby conductors, etc. In this section we will point out how the presence of this sort of noise could be indicated 
by lifetime measurements, and how it can be distinguished from the other sources of noise we have been considering 
in this paper. 

Electrical noise from a random field E (t) can produce spin relaxation if there is spin-orbit coupling present. 
Relaxation can occur by two distinct mechanisms: (1) spin-orbit-mediated virtual excitation to higher orbital states 
with spin flip and (2) modulation of the Rashba field. These correspond roughly to the Elliot- Yafet and D'yakonov- 
Per'el mechanisms in bulk. 
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For the first mechanism we have a Hamiltonian 

H = --g^ B B-a-eE(t).f+a[a x --a y -yp[a x --a v - 

This produces relaxation that is physically analogous to the spin-phonon mechanism, and the derivation is parallel, 
so we omit it. We obtain 

^ 2 / /|2fc4 



_L_ = ^em t g^ B B) ^ ^ bB/K) (^^j (a + pf (l + cos 2 9) (1 - sin 2$) a S Em (gn B B/h) B 2 L 8 



where Se x (gf^BB/h) is the spectral density of the E x autocorrelation function, evaluated at the qubit operating 
frequency, and L is a measure of the diameter of the dot. r] is a numerical factor of order one that depends on the 
shape of the dot. We write the applied field B as B = B (sin cos <3>, sin sin <£, cos 0) . = is the z-axis of the lab 
frame. Note that, as before, the relaxation rate decreases as the excited state splitting increases since the spin mixing 
which allows the electric field to relax the qubit is via the excited orbital state. 

The second mechanism is physically distinct in that it does not involve orbit ally excited states; instead the noise 
is converted to random time-dependent effective magnetic field on the spin. It is sufficient to consider a Rashba 
Hamiltonian, 



1 f d d 

-g^i B B - a - ect\E z (t) la x — - a y — 



which leads to 



Y = v'tftf ( a o) S Ez (g^ B B/h) , 

where again rj f is a geometry-dependent constant of order unity. Here, the rate increases with smaller dot sizes 
(L 2 = ^^). Reasonable values for the parameters are = 10 -5 nm 2 [26J , and L = 50 nm; evaluating this formula 
leads to 

^ ~ 0.15" 1 x S Ez (g^BB/h) , 

if Se z (qiabB) is measured in V 2 -s/m 2 . Zimmerman et al. determined the strength of electrical noise in a Si SET 
structure by measuring fluctuations in the peak separations of Coulomb blockade oscillations, but so far this type of 
measurement has been peformed only at frequencies much less than 1 GHz, which makes it difficult to estimate the 
noise magnitude at typical qubit operating frequencies in real structures. But the only Independence in T\ comes 
from Se z (9/abB/H) , which is likely to vary extremely slowly with B for any mechanism that one can think of. This 
means that defect-dominated electrical noise can be easily distinguished from other relaxation mechanisms by the fact 
that it is ^-independent. 

Magnetic noise from quantum and thermal current fluctuations in metallic portions of the circuit will produce a 
fluctuating magnetic field at the qubit that can relax the spin. No spin-orbit coupling is required for this mechanism 
to operate. This effect has recently been calculated by Langsjoen et al. [54]. These authors found values of T\ 
of order seconds for typical quantum dot architectures. The field and temperature dependence is given by \jT\ ~ 
B coth {ii B B /2k bT) , which reflects the photon density of states and the Bose function. The field and temperature 
dependences are again distinctive. 



D. Spin relaxation due to nuclei 



Hyperfine coupling of the electron spin to nuclei can give a very small admixture of the opposite spin state into a 
predominantly up or down state. This mechanism would give a T\ that depends relatively weakly on field. However, 
theoretical estimates give a small magnitude for this effect [12j[T9j[55]. This conclusion would of course be strengthened 
in isotopically purified Si 28 . Furthermore, this mechanism is not specific to dots and should occur also in donor spin 
relaxation, where T x w 0.25 x 10 4 s for Si:P at B = 0.32 T and T = 1.25K [U [12]. It does not seem to have been 
observed. Hence this mechanism is probably negligible at the fields and temperatures under consideration here. 
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E. Two-phonon processes 



At higher temperatures there is an activated two-phonon contribution from SOC mixing in silicon quantum dots. 
In the Si:P system, these Orbach processes dominate for T > 2K [56]. We can use the methods proposed by Castner 
[14J to estimate the Rashba + Orbach spin relaxation path. We arrive at 

M 2 so Y^ g n(E lg ) 



rbach 



where we estimate for a circular dot that Mso ~ (1 1| x \d 1} I (1| x \d) ~ 2 (a + f3)^/m* / Ei g as the spin-mixing of the 
two states and Y is the orbital relaxation time of the first excited state. For a circular dot with a parabolic potential 

r xy K Z exp [- 1 A| /kT] . (15) 



T Orbach *V h 4 7Tpsi 

Thus the temperature dependence of Ti at higher temperatures can give an accurate measure of A, a technique already 
used to find energy splittings of donor states. This can provide a check on transport spectroscopy determinations of 
this quantity. 



VII. VALLEY RELAXATION 



Electrons in lateral silicon quantum dots typically reside in the two degenerate conduction band minima along 
the z direction. This doubles the number of levels in the dot relative to the r-point-centered, direct band-gap III-V 
quantum dots as was described in Section [TTJ Here we wish to consider the relaxation times of these excited valley 
states as we have done for low-lying orbital and spin states above. Castner was the first to calculate the relaxation 
across different valley states from the 2p to Is levels in donors [13J. This has been repeated in Ref. [4 J for Li donors 
and for P and Li in strained silicon in Ref. [30]. A similar calculation can be done for lateral quantum dots where the 
interface in z is assumed perfectly flat and smooth, and thus the valleys can be considered good quantum numbers 
in the usual Kohn-Luttinger approximation (and the problem is separable in the three dimensions). We first consider 
this ideal (or "ID") case (which may be relevant in some experimental situations) and then comment on the more 
usual case of significant valley-orbital wave function mixing due to imperfect interfaces. 



A. Ideal (ID) interfaces 

We are concerned with relaxation across the same orbital and spin states but between valley states (v = +/— ) in 
a silicon quantum dot, particularly the relaxation of the lowest excited valley state with no change in spin or orbital 
number (type 3 in Figure 2|. O ur approach to valley relaxation follows the same procedure as exact orbital relaxation 
(Section |V| and Appendix IX A), where in this case we replace the matrix element with the inter- valley matrix element: 



M + _ = (ns+| H ep \ns-) . 

Assuming no valley-orbit mixing with higher states (separable wave functions in (x, y) and z), the wave function for 
an electron in a lateral silicon quantum dot reads 



^)(r) = F x , y (x,y)F x (z) a?>u z (r)e^ z + 



(v) 

ar'u- 



5 (r)e- 



-ikry 



F X) y(x,y)F z (z) 



G z 



C Gz exp (i(G z + k m )z) + a 



K- z 



C K _ Z exp (i(K_ z - k m )z) 



where k m is the location of the minima along the z-axis, a^(z^-z) = l/<\/2{l, ±1} (though these may be 
complex in the general case), and we have expanded the Bloch function in reciprocal lattice vectors, G or K, 
^( r ) = ^G^G ex P[^G • r]. The first five terms of the Bloch expansion contribute 90% of the wave function am- 
plitude (values from a recent study are listed in Table 1 of Ref. |57J). Here, the envelope functions of the two states 
are the same, the spin states are the same, but the Kohn-Luttinger oscillations are out of phase (see Figure [8]). We 
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Figure 7: Longitudinal and transverse phonon valley relaxation rates for ideal interfaces versus valley splitting, A v . The dashed 
lines are with the height of the wave function, zo, held constant (note the peaks at 23 meV and 11 meV respectively) and the 
solid lines are the relaxation rates with the changing extent of the wave function in z (due to the electric field) included. The 
vertical dimension of (ideal) triangular quantum well wave function as a function of valley splitting (theory from Ref. |31|). 
Note that the tighter the confinement in z, the faster the relaxation (due to broadening of th e wav e function in momentum 
space.) These numbers are for Uki = 1 and should be renormalized down by u\ x (see Appendix IX B). 



assume that the wave function consists of Gaussians in all three dimensions. Following our exact orbital relaxation 
calculation, the valley relaxation rate of a parabolic, circular quantum dot in a [001]-strained silicon quantum well is 



x 2 exp(-z (A v )ql/4) 

V°/coJ _ 



exp 



-T 2 A 2 

4 h 2 v 2 



h 3 vf 



exp 



-r 2 A 2 

4 h 2 v 2 



" 2 (P 2 



(ZIP? 



(16) 



where 



and 



P™ = J x n exp [A s x 2 ] sinh 2 [B a x] dx 



A s 
B„ 



-2 r 2 
^QAs l x o 
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Figure 8: Results from ID tight-binding model for wave functions for ground and excited valley states and valley splitting for 
a 10 nm quantum well. The z-component of the electron dot wave function is the output of a 2-band tight-binding calculation 
(points) which has been interpolated (line) for a typical SiGe heterostructure with a quantum well of 10 nm, barriers of 150 meV, 
and a large growth direction electric field due to space-charge separation from the donor layer of 6 x 10 6 V/m. Valley splitting 
in realistic silicon quantum dots will likely be reduced versus the ID results presented here due to interface roughness/steps, 
etc. 



where A v is the valley splitting, zo is the extent of the wave function (assumed gaussian) in z and q u is the phonon 
wave length of the emitted Umklapp phonon, q u = 0.3k max . The details of this calculation are given in Appendix 

At first glance, the valley relaxation rate has a A 3 



Let us compare Eq. 



16] to pure orbital relaxation, Eq. [10 
dependence as opposed to a A 4 in the orbital case (assuming parabolic dot potentials for both and matrix elements 
given due to gaussian wave functions). To understand this remember that for valley relaxation this transition occurs 
within the lowest manifold (both initial and final states have the same s-like envelope function) such that the matrix 
elements M ~ 1. In the orbital case, we must calculate matrix elements from 2p-like to ls-like states, such that 
M oc Xq oc A. The valley relaxation expression also includes prominently a exp(— z^q^/A) prefactor absent in the 
exact orbital case (Eq. 10). This prefactor predicts that the phonon relaxation rate will be peaked at the Umklapp 
phonon energy (assuming zq is constant with A v , which it isn't). Equation [l6| also shows the importance of the zq 
extent of the wave function; decreasing zq increases the relaxation rate. These effects are related, in that Umklapp 
phonons at q u = A v /hvij which connect valleys in neighboring Brillouin zones are the most efficient relaxation channel 
(see Appendix IX B for more details) . Figure [7] explicitly shows the valley relaxation rate in the two cases of fixed zq 
wave function height and wave function height that changes accurately with electric field and valley splitting. It turns 
out that the Bloch coefficients to the nearest valley at 0.3k max are most efficient and phonons are then emitted in the 



21 




O K 



realspai 



— o 

o> — O 

k - space 



Figure 9: Illustr ation of Umklapp phonon process which enables valley relaxation in ideal, silicon quantum dots (see text and 
AppendixjlXB). 
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Figure 10: Total valley relaxation time (due to both transverse and longitudinal phonons) for an ideal-interface, quantum dot 
excited valley state in silicon versus valley splitting, A v . The exact results including all multipole contributions are given by 
the blue solid lines while electric dipole approximation (see text) results are given by the purple dashed lines. Typical valley 
splittings in silicon quantum dots are less than 1 meV. 
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z direction. As zq gets compressed, not only does the valley splitting increase due to interface scattering (approaching 
the "critical" Umklapp phonon energies at 13.4 meV (longitudinal) and 23.2 meV (transverse)), but the wave function 
gets broadened in momentum space, allowing lower energy phonons to connect the two opposite valleys (see Figure 
[9]). So Umklapp valley relaxation is possible even at valley splittings smaller than q u . This effect causes the relaxation 
rate to continually increase as the valley splitting approaches q u , while the other "typical" exponentials kick in at higher 
splittings to cause a bottleneck effect as in the orbital case. The line width should be weakly dependent on the size 
of the dot in the lateral dimensions (as is the valley splitting) and much more so dependent on changes in z extent of 
the wave function. 

The valley splitting varies roughly linearly with E-field in a perfect quantum well where the electron only sees one 
side of the quantum well. Figure [7] shows the valley relaxation for an ideal interface as a function of valley splitting 
with a zq that changes correctly with A v . We account for the change of wave function size in z as a function of valley 
splitting. We take the ideal theoretical maximum valley splitting as (given in J) 

A„ = ^ « 2.3 x l(T-£ 

where E (V/m) is the electric field in the z-direction, AE C is the conduction band offset, v v — 7.2 x 10 _11 A£? C (in 
eV m with AE C in eV) [31J. Thus, the valley splitting A v depends on the i^-field, which also determines the extent 
of the wave function in z. zo now is a function of the E-field in z (which varies by device and can often be changed 
somewhat in a single device). For this we define the wave function in z as (assuming a triangular potential): 

^ = 1.4261v^ A%[kz - 2.3381) 



where k = l/z = y / 2m*eE/h 2 (in meters). Now we replace E with E(A V ) = A v /(2e x 7.2 x 10 n ) (in V/m). Ai is 
the Airy function. So the extent of the wave function in, zq(A v ), changes with the valley splitting as: 

, A , /ft 2 (7.2xl0- n )^ 1/3 



m*A v 

For completeness, we may also look for an expression for the valley relaxation in the electric-dipole approximation. 



A reasonable approximation is to set q\ s = in Equation [16] which, to leading order for a parabolic potential in all 
three dimensions, gives 

^"^■^(-iA.CA.)')^^. (17) 



The exact and approximate valley relaxation times are compared in Figure 10 



B. Comments on imperfect interfaces (tilt, roughness, and alloy composition) 

Our expression, Eq. [l6j for valley relaxation is for perfectly smooth interfaces. As we discussed above, imperfect 
interfaces will cause mixing between the orbital and valley states. In realistic devices, miscut, alloy variability, surface 
roughness, etc. will be present. This distorts or mixes the orbital and valley states such that valley is no longer a 
good quantum number [31] (e.g., the envelope functions within the s-manifold states can now be different and/or 
non-5-like). In this case, one must realize that the wave functions of orbital states will be different from the Gaussian 
wave functions assumed for s-like and p-like dot states used in the orbital relaxation section. In reality they will be not 
be separable, and there will be sample dependence. However, all these energy relaxation calculations are proportional 
to M 2 , where M is the dipole matrix element. In contrast, the energy gap dependence is much higher: A 5 (A 4 for 
parabolic dots) and the field dependence is B 7 . A is much easier to determine than M and is much more important, 
which means that our predictions are still useful. A likely exception is our ideal calculation for possible long-lived 
valley states. In the non-ideal case, there is likely to be valley-orbit mixing, providing an avenue for relaxation via 
G = phonons as in the orbital relaxation case. While a full theory of this mixing is possible (utilizing appropriate 
wave functions, e.g., following [3Tj l33j), the matrix elements depend on the exact specification of the interface for the 
dot being measured (which is difficult to ascertain) and are not considered here. Because the orbital relaxation is 
fast, the long-lived predictions for valley excited states will be wrong in this case. Although the degree of the this 
mixing/distortion depends on the specific device in question, we can crudely write that 

perfect = ^.p^ + (J _ T ideal 

where f m i x is the notional fraction of orbital wave function mixed in with the valley wave function and Tio (A v ) is 
the orbital relaxation rate across the measured valley splitting of the measured state. In some cases f m i x could be 
calculated [31 J but generally it will depend on microscopic details of the silicon quantum dot. 
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VIII. COMPARISON WITH EXPERIMENTS AND CONCLUSIONS 



Some spin relaxation experiments have been performed on electrons in silicon dot devices (see Figure [5|. Our 
theory for spin relaxation in the electric-dipole approximation should be valid only well before any degeneracy between 
Zeeman and orbital splitting is reached, the regime of a silicon spin qubit. In Ref [47J, T\ was measured in a lateral, 



depletion mode Si/SiGe quantum dot for three different magnetic fields between 1 and 2 T. Comparing to Eq. 12, the 
data appears to follow a B~ 7 rate but with relatively few points it is difficult to make firm conclusions. These data 
points fall on the 0.4 meV orbital splitting line (obtained using the value of SOC from Ref. [46J), which is roughly 
consistent with the lithographic size of the dot but differs from the A = 2 meV value used for the theory model in 
that paper (though it is unclear if the experimental orbital splitting was measured) [47]. T\ has also been measured 
in a metal-oxide-semiconductor system with the dot near the Si/Si02 interface [48J; they report measuring an orbital 
splitting in this device of roughly 0.4 meV. There are five data points that fit a Ti ~ B~ 7 law reasonably well; however, 
this behavior is observed only for B > 3 T. Here theory predicts a roughly order of magnitude shorter spin relaxation 
time for a dot with A = 0.4 meV (the data > 3 T behave as if the orbital splitting is 1 meV). For B < 3 T, Ti is 
roughly independent of B, with T\ = 40 ms (possibly limited by charge noise or some other mechanism). Finally, a 
recent measurement of T\ in a laterally gated Si/SiGe dot in a doped device yielded T\ = 2.3 s at a field of B — 1.85 
T. None of the data shows evidence of hot spot (fast spin relaxation cusp) behavior. 

On the whole, these experiments give good evidence that the SOC-mediated spin-phonon interaction is the dominant 
channel at high fields: both the B dependance and, importantly, the overall magnitude are consistent with theory. The 
value of the Rashba coefficient that is used in Figure [5] may not be appropriate for a MOS structure (although inversion 
layer spin relaxation measurements show relaxation times [58] within an order of magnitude of those found in SiGe 
QWs) and can vary from device to device based on material, interface roughness, electric field at the interface, etc. 
[26| l5T]. A much smaller SOC constant might explain the order of magnitude difference between theory and experiment 
for the various systems. It would be very useful in the future to attempt to characterize the SOC strength in these 
wafers by other means, for example via 2DEG spin relaxation |20| . although the SOC may vary at a microscopic level. 
It is also very important to check the dependance of T\ on field direction. This has not yet been done experimentally. 
Looking beyond the electric-dipole approximation, there may be a strong and sudden increase in the relaxation rate 
for spin relaxation when the spin splitting matches the orbital splitting; no such effect is expected if the 1st excited 
state is a (pure) valley state. A complication to this picture may be a lack of mixing with the orbital state if the 
SOC is only Dresselhaus-like [52J (even in the ideal case). Lastly, our assumption of Fock-Darwin wave functions may 
not be correct (do to dot asymmetry, interface roughness, etc.); leading to matrix elements between dot states that 
influence the relaxation rate up or down. 

While these results are encouraging, they do not constitute a complete vindication of theory. Hence we discuss 
how to combine the results of different measurements to fix some universal quantities. Specializing Eq. [6] to the first 
excited state, neglecting anisotropy and the dipole moment in the z-direction, noting that v^ 7 >> vj 1 ', using Eq. [TJ 
and taking ksT « A, we find an orbital relaxation rate from the first excited state to the ground state: 



r 12 = 



4IAI 5 - 2 , 



Mi 12 ) 



2 



(18) 



105 h 6 7rp S ivJ 

With similar assumptions for the spin relaxation and specializing to an in-plane field, we have from Eq. [12] 

(12) 

We wish to eliminate the poorly determined quantities Mi ; and £, both related to the size of the dot, and the 
electron-phonon coupling strength and phonon velocity, in favor of measurable numbers (as far as possible). To do 
so, we use Eqs. [8j and[T3] 

4 A 8 h 2 

Here Ez = 9I^bB is the Zeeman splitting. A is the energy of the first excited state, which can be measured 
independently, by transport spectroscopy or looking at Orbach processes at higher temperatures. /3 can also be 
measured by other means, though this is not straightforward [59J. Thus it would take a combination of measurements 
to use the absolute magnitude of T\ as a test of theory. Note that this analysis assumes that the first excited 
state is purely orbital in nature and that the fast orbital relaxation time can be measured. The former may be 
ameliorated in the valley case where there is strong valley-orbital mixing, or if not, Y21 can be replaced with the 
ideal valley relaxation time. The latter may be a difficult experimental constraint. Failing this, the 1/Ti ~ B 7 
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1/Ti ~ [1 -f 2n# (Ez/kBT)} , and the field anisotropy given by Eq. [l2| provide the best tests. We note also that 
the field anisotropy of the spin relaxation, Eq. [l2j can help determine the relative contributions of the Rashba and 
Dresselhaus-like SOC contributions. 

We have focused on the lowest lying state configurations of a silicon quantum dot that are most relevant for quantum 
computing and have found that, in general, spin-based quantum computing benefits from orbital and valley states 
being as high in energy as possible. We began by considering phonon relaxation of excited orbital states across the 
same valley state in a lateral silicon QD. We found that orbital relaxation could be dramatically faster in biaxially 
strained silicon than in the bulk. This, for example, speeds up spin qubit initialization via optical pumping schemes 
[27\ as well as possible leakage to excited states via phonon excitation. The phonon bottleneck effect will eventually 
decrease the orbital relaxation rate but only for unrealistically small dots. In contrast, spin relaxation can easily be 
seconds (even when the magnetic field points along the growth direction) and T\ increases for small dots and low 
magnetic fields. At small magnetic fields, charge noise could play a dominant role. Valley relaxation can also be long 
in ideal dots, especially for small valley splittings. In non-ideal dots, although not quantitatively considered here, 
valley relaxation will likely be comparable to orbital relaxation. 

The theory proposed here depends on the correct identification of excited states as either orbital excited states or 
valley excited states. Theoretical considerations for spin and valley relaxation in cases where the states are not purely 
orbital, valley, or spin — in other words they are mixed due to, for example, disorder or surface roughness in realistic 
devices — and for regimes beyond small B-field (where degeneracies come into play) are subjects for future work. 
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IX. APPENDIX 



A. Orbital relaxation: beyond the electric dipole approximation 

To calculate the matrix element of H ep between orbital states m and n in Eq. [5] beyond the electric dipole 
approximation, we begin with the full expression (see Eq. [4|: 

(n\ H ep \m) = / ( H ep) ^ndV 



(21) 



We proceed following the derivation by Castner [14] . This derivation will be useful when valley relaxation is considered. 
A function which is periodic with the period of the lattice may be expanded in a Fourier series in the reciprocal lattice 
vectors Q^, so 



and the integral in M mn becomes 



FA 2 e-^'+^+v-^dr. 



The envelope probability can be Fourier transformed, 

IW! a =(2^E/0Oe*"- r . 



Plugging this into the integral in M mn gives 
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which equals 

£ Cfc-fc, (iE/ ( k ") 53 ( k > - k,' - q + Q, + k") = £ Cki-k^i-* + k, + q - Q„). 
Finally, the matrix element is given by 

(n| ff ep \m) = <A [3d (e s • q) + 3„ (q • K«) (e s • K«)] ^ £ C£_ fe . (-k, + k, + q - Q„). (22) 

We are calculating an intra-valley scattering process (orbital relaxation with no change in valley state) so Q u = 
is the dominant term, lq = kj , and a m — a n which gives 

(n\ H ep \m) = cL< [S«i (e (q, A) • q) + S u (q • K«) (e (q, A) • K«)] Cj^ / m "(q), 

and for the most relevant transition, 

(2| H ep |1> = £<,< + S u e z(Zz ] C£_ fe . / 12 (q) 
= [E dqi +E u e z q z }Cl_ k . / 12 (q), 

where qi = qA = A/hvi for longitudinal phonons but qg = for transverse phonons (</a = A/hvt). C^._ k . is the first 
coefficient in the Bloch wave expansion (see Table X for the largest contributions). 

The envelope function of the ground state QD wave function in the absence of a magnetic field in the lowest approx- 
imation is a product of Gaussians, F^\r) = F(x,y,z) = F(x)F(y)F(z), where F(x) = (2/tt) 1 ^ 4 x^ exp (— x 2 /xq) , 
F(y) = (2/7r) 1/4 2/o" 1/2 exp (-y 2 /yl) , and F(z) = (2/tt) 1/4 Zq 1/2 exp (-z 2 /zl) . The excited state, if y > x » z , 
is F^(r) = F(x,y,z) = F(x)F {2 \y)F(z), where 



F< 2 >(j,) = [2/^/yij (2/^Vxp (V/^) • 
Then, the overlap integral is given by 

/( 12 )(q) = J F( 2 )(r)e <q ' r F( 1 )(r)dr, 
/ (1) (<fe)/ (1 %y)/ (1) fe) = y F^ixfe^dx f F { - 2 \y)e i ^F^{y)dy J F^(z) 2 e^ z dz 

( 1 2 2\, %) M 2 2^ / 1 2 2^ 

= ex P ( -g^ofe J ^-9yexp I --J/offy ) ex P ( -g^0<7 z J • 



Inserting these into the Golden Rule, we find that 
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r 2 i = ^]T|M 21 | 2 o(A-fia, q>s ) (23) 
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continuing, 
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We are left with calculating the three angular integrals / which have units kg 2 /s 2 and are defined as 
L = [ \[E d q Ai ^E u e z q As q z ] f 12 (q As ci)\ 2 d6d<l) 

J V s 

= J sin6> ^ [-dQAi + Z u e z q As q z ] 2 exp (^-^xlqi^j ^^exp ^-^VoQ^j exp (-^o<? 2 
We can immediately point out that I t2 = because e z (t\) = 0. Then, 

L = J sin ^2° [^dQAl + Z u e zs q As cos 0} 2 exp ^-^o^L sin 2 (9 cos 2 
x sin 2 #sin 2 </>exp ^— -^/q^As sm2 # sm2 <^ ex P ^ — ^o^As cos2 ^ dOdcj). 
If we assume an approximately circular dot in x and then using cos 2 + sin 2 <fi — 1 we can do the integral easily, 

.2 af 2 /-I 2 \ 

As^ ) 



t ( ^ 2 2 \ f QAsVo r— i " l 2 /i 2\ (^(2 2 

I s = Trexp (^--x (7 As J y [XdQAi + ^e^AszJ (1 - x ) exp I -(aj - z 

It = -exp (-^§«L) «M [~ 2 ,L(A? - A 2 ) + 2~,~ n ,L(A 2 - A? ) + 5^(Af - A? 
^ = ^xp (-\x 2 q 2 A ^ %#S 2 ,i £ [A 2 - 2A\ + A?] , 



where 



Finally, the orbital relaxation rate for a parabolic dot (in all three dimensions) from its first excited state is given by 
(with E21 = A, the common notation) 

21 2(27T)2pSiJi 4 

I exp(-^ g L) A*_ [Hg( ^ _ ^ + 2Ed ^ {Af _ Af) + ^ _ Af)] (25) 

[ V l a V l 

+ ex P (-y g L) ^ s , [A? _ 2At +A ,y (26) 

This reduces exactly to the expression for orbital relaxation within the electric dipole approximation, given in Eq. [6j 
when q\ s = 0. 

B. Valley relaxation (ideal case) 

1. Valley relaxation in a three-dimensional parabolic quantum dot 

We consider valley relaxation in a lateral silicon quantum dot. We begin where we left off in our exact consideration 
of orbital relaxation. Our expression for the electron-phonon matrix element, Equation [221 was 



run 



<< [3d (e (q, A) • q) + E u (q • K«) (e (q, A) • K<*>)] a* £ C^/'^-k* + k j+ q- Q„), 



where / is the Fourier transform of the envelope function overlap integral. Since the valley transition involves a change 
in crystal momentum, there are no intravalley terms from this expression and we must consider high wavenumber 
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phonons that can connect the two valleys which are separated in the first Brillouin zone by 2k m i n . In this case, i j 
and kj = — and the matrix element becomes 

M mn (i ^j) = J2 a m°i [ S <** + 5„ (q • K«) (e (q, A) • K«)] a* £ (-2^ + q - Q„), 



where = <?a for longitudinal phonons but zero for transverse phonons (q A = A/hv s ). The matrix element will only 
be large for values of q ~ 2k^ — K^. The shortest wavenumber phonon to connect the two valleys is the Umklapp 
phonon across the Brillouin zone where \q u \ = 2(q rnax — qmin), where q m ax — ir/asi- Thus, for the +z and —z valleys, 



M as = - [E dQl + E u e z q z ] a* q Cl Q [f {1 \-q u z + q) - f W (q u z + q) 



(27) 



This is just the result of Castner as a component of his calculation of Raman spin transitions in donors. The major 
difference between the donor and QD calculations (in the ideal case) are due to the different envelope functions 
(impurity vs. parabolic). We next require the Fourier transform of the QD envelope function: 

/ (1) (q) = f F^(r) 2 e icL - r dr 

f il) (qx)f il) (qy)f il) (q z ) = J F^\x) 2 e iqxX dx j F^\y) 2 e iq y y dy J F^(z) 2 e iQzZ dz 



exp 



8 



2 2 

x q x 



exp 



8 



vWy J ex P 



- 2 2 

; z oQ z 



Again we have considered the case where the z dimension of the wave function can be approximated as a simple 
Gaussian (which is for our consideration a good approximation). 

Looking at Eq. 27, we see that the /-functions are heavily peaked at q u 0.3q max = 0.3ir/a — 1.74 • 10 9 m _1 ) 
in the z direction and at in the x and y directions. Since phonons of this magnitude are needed to connect the 
two valleys, resonant phonons close to this will increase the matrix element leading to increased relaxation. However, 
slightly off-resonant phonons can also cause a transition due to the widths of the /-functions which broaden as zq 
gets smaller (see Figure [9]). In donors, the valley splitting tends to be around 11 meV, not far off of this wave vector. 
In silicon quantum dots, the theoretically predicted values of the valley splitting range from to 3 meV depending on 
the extent of the z wave function. 

To calculate the valley transition rate we employ the Golden Rule, 



2tt 
"ft 



sin 



\M a 



r<5 U 



A 

hv. 



(28) 
(29) 



where we have summed over phonons and the emitted phonon has wave number q A = A/ hv s . At cryogenic temperature 
there are absolutely no large wave number phonons, so we need only consider spontaneous emission. Incorporating 
our expression for the matrix element, we find (with F(q u , q) = [f^(—q u z + q) — f^\q u ^ + q)]) that 



2tt V 

K + i) {cp 

8(2ir) 2 p S ih 



q dq sin OdOdcj) 



1 



i [E d q t + E u e z q z ] a* q Cl Q F(q u , q) 



S(q- Qa) 



poo r ^ 

/ q sin 6^ \[E d qi + E u e z q z ] F(q u ,<i)\ 2 5 (q - q A ) d6d(j)dq 

e J0 J v s 



k ) 



Taking the delta function, the rate becomes 

(n„ + I) iut. x2 



n, the rati 

S(27r) 2 p S ih 



K + i) {ci 



k J 



S(2n) 2 ps z h 



^QAs / si 

^QAsL 



sinO— \ [E d q Al + E u e z q As q z ] F(q u ,q As q)\' 
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where q = ej. We are again left with calculating three angular integrals / which have units (kg 2 /s 2 ) and are defined 

as 



L=lMM = / / sin6>^ \Z d q A i + 'E u e zs q As q z \ (-q u z + g As q) - (<2u^ + #A s q 
./o ./o v s L 

it 2 =0. If we assume that the dot is circular , then the /-functions simplify, 



dOdcj). 



f {1 \±q u z + g A q) = exp ^-^o (<2a^) 2 - (<?a^) 2 ^ exp ^~^o (qaQz ± g n ) 2 ^ 
= exp \-\xlq\ sin 2 ] exp \-\zl (q A cosO± q u f 



Replacing q with its components we can further simplify /: 

r>2ir p7T 



r Z 7T n 7T j 

^=Mi,t 2 = / / sin# — 

^0 ./0 v s 



[^dQAi + ^ u e z q As cos 0] exp --a g A sin 6> {} 



d0# 



where 



= ^ y sin #^2 + ^uQAsZz cosO} 2 exp ^-^o^a sin2 ^ {} 2 dOd<l) 

{} 2 = jexp ^-^o (QaicosO - q u f^j - exp ^~^o (qai cos + g n ) 2 ^ | 
= exp ^-i^g 2 ^ exp ^~^o [<?icos 2 0]^ 4sinh 2 Q^2<? A g n cos 0^) . 



So, taking the trivial integral (no e zs depends on 0), 

I a =iMte = 87rexp (~^*o^ y sin0d0 ^2 1-dQAi + S n g As e 2 cos6>] 2 

x exp ^— ^o^az sm2 ^ ex P ^ — [#A cos2 sinh 2 ^ZQ2q/\q u cos O^j . 

Now, we explicitly consider the integrals for s = I and s = £2- Evaluating the longitudinal and transverse (£2) integrals 
separately: 

7/ = 87rexp ^-i^ 2 ^ exp (-^x 2 q 2 Al ^J ^ 

x y <£z [E 2 d + 2~ d ~ n x 2 + S 2 x 4 ] exp QtfL [x 2 , - ^] x 2 ^j sinh 2 Q^2<? A ^ n ^ 



and, similarly, 



1 2 2^ / 1 2 2 \ # A t 



7 t2 = 87rexp ^--z o9u j exp ^--a^At J — ! 

X y ^~ 2 - ^] ex P Q^At [^0 - *o] Sinh2 Q^0 2 <?At<2u 

The integrals have no analytical solution so we define a numerically tractable integral function 

P s n = y x n exp [Ax 2 ] sinh 2 [Bar] 
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where 



A = 4<?L [ x l - z o] 



and rewrite our solutions: 



Ii = 87rexp ( --^o^n ) ex P 



1 



7 x oQAl 



2 2 \ QAl p i ^ p2 I ^2 p 41 



and 



i"t 2 = 87rexp ( ~-z q u ) exp 



1 



» x oQAt 



2 2 \ QAt^2 \2 p 2 _ 4 p 4 



Finally, the valley relaxation rate of a parabolic, circular quantum dot in a [001]-strained silicon quantum well is 
(n q = and C = 1) 



Tas = (n q + 1) {Cl) [ 



exp 



(=4*) 



exp 



2 2 \ 3 

X Q^A^ \ /~2 p i ^ p2 I ^2 p 4\ 

A J ~^2 l^dM + ^d^u^l + J 



r>2^2 
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